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o Abstract 

In the present work, we analyze the evolution of the scalar and tensorial perturbations and the 
quantities relevant for the physical description of the Universe, as the density contrast of the scalar 
perturbations and the gravitational waves energy density during the Bose-Einstein condensation 
of dark matter. The behavior of these parameters during the Bose-Einstein phase transition of 
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dark matter is analyzed in details. To study the cosmological dynamics and evolution of scalar 
+2 and tensorial perturbations in a Universe with and without cosmological constant we use both 

analytical and numerical methods. The Bose-Einstein phase transition modifies the evolution of 
i— i gravitational waves of cosmological origin, as well as the process of large-scale structure formation. 

^ 1 Introduction 
O 

In what concerns the standard cosmological cenarium (ACDM model) the observational data [U [2] 
show that the ordinary baryonic matter is responsible for approximately only 5 % (neglecting radiation) 
of the total energy density of the Universe. Nearly 70% of the total density is the contribution of the 
so called dark energy (DE), responsible for the observed acceleration of the cosmic expansion [31 0]. 
I The other 25% is in the form of a weakly interacting matter, called cold dark matter (CDM) or simply 

dark matter (DM), whose properties are inferred from its gravitational interaction with visible matter, 
k> radiation and its effect on the large scale structure of the Universe. 

The existence of DM was initially observed due to discrepancies between the calculated mass of the 
observed matter and the mass determined from the gravitational effects. Subsequently the observation 
of the rotational curves of galaxies, the velocity dispersions of galaxies, the gravitational lenses and 
the large scale structures [51 [6] suggested the existence of the DM. There are several hypothetical 
candidates for DM particles e.g. weakly interacting massive particles (WIMPs), axions, hidden dark 
matter, sterile neutrinos and light gravitinos [7J and although the nature of DM is still unkown a great 
effort has been made in an attempt to detect the DM particles EJ [7J . 

The BEC process, that is a very well observed phenomenon in terrestrial experiments, occurs when 
a gas of bosons is cooled at very low temperatures, near absolute zero, what makes a large fraction 
of the particles occupy the same ground state. This kind of condensation could have occured at 
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some moment during the cosmic history of the Universe. The cosmic BEC mechanism was broadly 
discussed in [51 [9]. In general the BEC takes place when the gas temperature is below the critical 
temperature T crt < 2irh 2 n 2 / 3 /niks, where n is the particles density, m is the particle mass and Ub 
is the Boltzmann's contant. Since in an adiabatic process a matter dominated Universe behaves 
as p oc T 3 / 2 the cosmic dynamics has the same temperature dependence. Hence we will have the 
critical temperature at present T crt = 0.0027 K if the boson temperature was equal to the radiation 
temperature at the red-shift z = 1000 [9|. During the cosmic adiabatic evolution the ratio of the 
photon temperature and the matter temperature evolves as T r /T m oc a, where a is the scale factor 
of the Universe. Using as value for the present energy density of the Universe p = 9.44 x 10~ 30 g/cm 3 
BEC will happens if the boson mass satisfies m < 1.87eV. 

Despite the great sucess and simplicity of the standard ACDM cosmological model, it still remains 
as a phenomenological model and we can not explain why the cosmological constant density has a value 
close to the critical density today and why the measured cosmological constant is smaller than the 
quantum vacuum energy density by a factor of 10~ 120 . At galactic scales, this model has difficulties in 
explaining the observed distribution of DM around the luminous one. The simulations in the ACDM 
model predict that bound halos arround galaxies must have very characteristic density profiles with a 
well pronounced central cusp [10] , but the observations of rotation curves show that the distribution 
of DM is shallower than the predicted and it has a approximately constant density core [llj. The so 
called core-cusp problem can be solved by presuming that DM particles are ultralight scalar particles 
with masses of the order of 10~ 22 eV initially in a cold BEC [12J. 

Assuming the hypothesis that the CDM in a galaxy is in form of BEC the density distribution of 
the static gravitationally bounded single component BEC DM is given by p(r) = p* smkr/kr, where 
p* = p(0) is the density in the center of the condensate and k is a constant. Giving the conditions 
p(R) = and kR = tt, where R is the condensate radius, the condensate DM halo radius can be fixed 
as R = Trffils/Gm 3 ) 1 / 2 [13] where l s is the particle scattering length and m is the the particle mass. 
The calculated total mass of the condensate DM halo is M = A-K 2 (h 2 l s /Gm?f/ 2 p* = 4R 3 p*/ir [IB] . 
So the mass of the particles of the condensate are [13] 

™=(w) ss6 - 73xl ° (rinj (nsO eV ■ (1) 

Recently the cosmological BEC process of DM was investigated and in this model [H] it is assumed 
that the condensation process is a phase transition that occurs at some time during the history of the 
Universe. In this model the normal bosonic DM cools below the critical condensation temperature, 
that turns to be favorable to form the condensate in which all particles occupy the same ground state. 
In this new period the two phases coexist for some time until all ordinary DM is converted into the 
condensated form, when the transition ends. The time evolution of cosmological parameters as the 
energy density, temperature and scale factor is changed during the phase transition process [14] . 

In this letter we aim to investigate how the generation of cosmological perturbations produced 
by the space-time expansion, and the parameter associated with them, would be changed during this 
Bose-Einstein phase transition of dark matter. The scalar perturbations are related to the formation 
of structures like galaxy clusters and galaxy superclusters in the Universe. The generally accepted 
picture of the structure formation in cosmology is the one where the Universe started off in an extremely 
homogeneous and isotropic state, with initial conditions provided by an era of accelerated expansion 
called inflation. The tiny primordial density fluctuations, generated during inflation from quantum 
fluctuations of the vacuum, would later grow under the influence of gravity and eventually collapse to 
form the structures that we observe today. In general, the evolution of the density of the different field 
components of the matter are calculated analytically or numerically depending on the complexity of 



2 



the studied model. These calculations are made in such a way that they may be compared with the 
observations of large-scale structure of the Universe like, e.g., Cosmic Microwave Background (CMB) 
anisotropy maps derived from WMAP data [H [2], Type la supernova (SNIa) surveys [3l 0] and the 
baryon acoustic oscillations (BAO) in the Sloan Digital Sky Survey (SDSS) luminous galaxy sample 

The gravitational waves are tensorial fluctuations in the metric of space-time. This particular 
perturbation is not explicitly coupled with the energy density and pressure of the matter of the Universe 
and does not contribute to the gravitational instability that generates the cosmological structures we 
see today. On the other hand this study is of great interest because it supplies the specific signature 
of the metric theory of gravity. These waves are predicted by Einstein's theory of General Relativity 
(GR) but they still have to be directly detected. Great efforts have been done in this sense and there 
is hope that a new generation of experiments in space may allow this detection [16] . The gravitational 
waves spectrum frequency extends over a wide range of interest, from 10 -18 -ff2i to W 8 Hz, depending 
on the sources that generate those waves. 

The present letter is organized as follows. The basic properties of the Bose-Einstein condensate 
dark matter (BEC DM) and the ordinary dark matter are reviewed in Section 2. The cosmological 
dynamics of the density perturbations is considered in Section 3. The evolution of the gravitational 
waves is derived in Section 4. We discuss and conclude our results in Section 5. 



2 Ordinary Dark Matter and Bose-Einstein condensated Dark Mat- 
ter 

Following the present data [IJ[2] we assume a flat homogeneous and isotropic Universe, whose geometry 
is described by the Friedmann-Robertson- Walker metric, given by 

ds 2 = c 2 dt 2 - a(t)dx 2 , (2) 

where a(t) is the scale factor of the Universe that describes the cosmic evolution, t is the cosmic time 
and c is the speed of light. The gravitational dynamics is given by the Einstein field equations 

R/J.U ~ l^dfJ-vR = c 4 . (3) 

We also consider the Universe filled by a perfect fluid, described by the energy-momentum tensor 

= (pc 2 + p)uV - pg^ , (4) 

where p is the density of the fluid, p is the pressure and g^ v is the metric tensor. 

Imagining the Universe filled by baryonic matter, radiation, dark matter and cosmological constant 
and neglecting any possible interaction between those components by assuming that the energy of each 
component is individually conserved, the equations of motion are 

a 2 8ttG . . , . 

"2 = —^-{Pb+Pr+Px + PA.) > ( 5 ) 

CI O 

a a 2 8ttG 
2- + ^ = — ^(Pr+Px+PA) , (6) 

Pi + 3(pi + ^)- = ■ (7) 



c 2 / a 



where H = a/a is the Hubble parameter and i = b, r, \i A means baryonic matter, radiation, DM 
and comological constant, respectively. 
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2.1 Non-condensated Dark Matter 



We will presume that in the early phases of the Universe the DM was made of bosons particles 
with mass m x and temperature T, originated in equilibrium and decoupled at a temperature Tp. 
Particles that decouple from the plasma in the primordial Universe at any temperature Tjj has its 
momenta redshifted such that the distribution function at any time after de decoupling is related to 
the value of the distribution function at the moment of the decoupling and it keeps an equilibrium 
shape in both non-relativistic and extreme- relativistic regimes |14| . In the non-relativistic regime the 
non-condensated DM pressure is [H] 

Px = <?° 2 p x . ( 8 ) 
where a 2 = (v 2 ) /3c 2 is the velocity dispersion and (v 2 ) is the mean square velocity of the particle. 
Appling the equation of state Q into the conservation equation ([7]) leads to 

P { x 

Px = (a/a )^ 2 ) ' (9) 

where p = p x ^ when a = do- Replacing the evolution equation (JoJ) for the non-condensated DM 
into the Friedmann's equation §5§ we find that the evolution of a Universe filled by baryonic matter, 
radiation, non-condensated DM and cosmological constant is given by 

a 2 2 / ft[ 0) ni 0) 4 0) \ 

= Hl \Vkf + + (a/ao)W) + ' (10) 

where Hq is the Hubble parameter and = S-^Gpf /3Hq is the density parameter, both at a = dQ. 

The density parameters satisfy the relation + fi^' + £l x °' + Q\ = 1 and ao can be considered the 
present day scale factor. Since the DM is non-relativistic the cosmic evolution is little dominated by 
the term a 2 . In the standard model a 2 = 0. 



2.2 Bose-Einstein condensated Dark Matter 

The BEC is formed when all particles in a Bose gas occupy the same quantum ground state. This 
happens at a very low critical temperature T crt ~ 2irh 2 p 2 / 3 /m 5//3 /cs [SO [2], where m is the mass of 
the particles, p is the density and is the Boltzmann's constant. The coherent state is reached when 
the temperature is low enough or the density is sufficiently high |14j . In this model we will assume 
that almost all the DM particles are in the condensate form and that the DM halos are made of BEC 
DM at absolute zero temperature. 

Only collisions between two particles at low energy are relevants in a diluate and cold gas. These 
collisions are characterized by the scattering length l s independently of the two-body potential |14j . 
hence the interaction potencial can be replaced by an effective interaction Vi(f ' — r) = X5(f ' — r) 
[T71 [T4"] where A = 4Trh 2 l s /m. The ground state features of BEC DM halos is described by the 
generalized Gross-Pitaevskii (GP) equation [T71 [J3] 

ih 94 *^ - ^V 2 <P(t,r)+mV(r)<P(t,r) + U \<P(t, r)| 2 0(t, r) , (11) 
at 2m 

where 4>(t, r) is the wave function of the condensate, m is the particles mass, V is tha gravitational 
potential that satisfies the Poisson's equation V 2 U(r) = 4irGp and Uq = \ixf?l s jra. We can use the 
Madelung representation of the wave function [8l El |T7j, which is 

^rVv^xe^" , (12) 
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where p(t, r) = \<j>(t, r)\ 2 is the density of the condensate and the function S(t, r) has the dimension of 
an action. This transformation will make the generalized GP equation (11) breaks into two equations 



~dt T x 



+ V-( Px v) = , (13) 



= -M£)-* v (£)- w « ■ (14) 

where Vq = -(h 2 /2m)\7 2 Jp^l p x is a quantum potential and v = VS/m is the velocity of the quantum 
fluid. The effective pressure of the condensate [8j H3] is given by 

and the dynamics of the fluid is determined by the scatering lenth Z s and the mass m of the DM 
particles. Using the conservation equation Q and the equation of state (15) we find 

Px + 3 Px C 1 + ^x) ^ = . ( 16 ) 

with w x = 2irh 2 l s /m 3 c 2 . The evolution of the BEC DM density is given by 

C 

Px= (a/ao) 3 -u x C ' 
where C is an integration constant. Ensuring that p x = p x at a = ao we will obtain 

Px = — 7 i P y^ > ( 18 ) 

uj x (a/a Q y - p 0x 

with 

.. J°) J°)o(°) 

^Xfa UJ xPcrt il X c lq \ 

where the critical density of the Universe today is p^) t = 9.24 x 10~ 27 Kg/m 3 . And the Friedmann's 
equation @ for the BEC DM is 

a 1 \ [a/ao) 6 (a/a ) 4 a; x (a/a r - p 0x 



where Q 0x = Po x J ' pSi- 



2.3 Cosmological dynamics during the condensation phase 

The best description of the BEC dynamics seems to be a first order phase transition [13], since it 
was demonstrated, by analyzing the temperature and density dependence of the chemical potential in 
several theoretical models that describe the thermodynamical transition from the normal phase to the 
BEC phase, that none of them predicts a second-order phase transition |14U18[ H5]. On the other hand 
a first order phase transition characterize a genuine mathematical singularity, but for the case of a 
ideal Bose gas confined in a cubic box it was shown that a system with a finite number of particles can 
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present a discontinous phase transition that characterize a genuine mathematical singularity, provided 
that the pressure is kept constant [T9] . 

During the BEC process the first thermodynamic condition that has to be satisfied is the continuity 
of the pressure at the transition point. So the equations of state ([8]) and ( 15 ) give the critical transition 
density [13] 

the critical temperature [U [T7] is 

_ 2irh 2 , ^2/3 . v 

CH ~ C(3/2)2/3 m 5/3fc B ) > l^J 

where £(3/2) is the Riemann zeta function. The scale factor of the Universe at the moment of the 
condensation is p3] 

2rf SS«<»>^ , (23) 



Q>crt 



a \c 2 a 2 m 3rcrt x , 
and the critical redshift is 

During the transition phase the pressure P = P cr ( and the temperature T = T cr i are constant and 
after the beginning of the phase transition the DM density decreases from p™ orm , when all DM is in 
non-condensated form, to p x when all DM is condensated, when the transition phase ends. So we 
introduce the volume fraction of matter in the BEC phase 



norm 



and we can write 



f(t) = H' rx , (25) 

» v ' n BEC _ n norm ' v ' 

r'X r X 



(1 + "*/(*)) , (26) 





norm 

~ Px 




n BEC 




Px 


P 



,norm 

X 



„norm 

• (27) 



When the condensation starts, at t = t cr t all the DM is in non-condensated form, i.e., p x (t cr t) 
■ x ' rrn and /(tcrt) = 0- At the end of the condensation all the DM is in BEC form, p x (tBEc) = Pv 1 



and f(tBEc) = 1- With help of the energy conservation equation Q we find that [T] 

a(t) = a crt (l + rf(t)y 1/3 , (28) 

r = ^ (2Q) 

l + Pcrt/p™ rm C 2 ' 1 ' 

where a(t cr t) = a cr t and the phase transition ends at 

l + z BEC = (l+r) 1 / 3 (l + z crt ) . (30) 
Using the Friedmann's equation ^ we find that the evolution of the volume fraction is given by 



1,1 ^ 1+r/ V/ (l + r/) + ^f^(l + r/)4/3 + 4 )(l + %/ ) + O A , (31) 



dr V r J V («cri/a ) 3 (a cr t/a y 
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where 4 0) =/ x , ,,,, 
the other hand P cr t/p 



/Pert and T 



norm 2 
X 



flo*- Since /O^ < p" orm we will have r £ (—1,0) and n x < 0. On 

n x [14J and 



<7 << 1 and from equation (28) we can approximate r 



if we neglect the radiation contribution to the energy density t 



re equation (31) has the solution 



/(*) 



rVL n 



1 + n c e~ 3( - QA / H o ^(t-tcrt) 
1_0 e -K^h/H^){t-t cr t) 



(32) 



where H 1 is the Hubble time and 



a. 



a 



(0) 



(a crt /a ) 3 x 



(33) 
(34) 



If we neglect both radiation and cosmological constant contribution to the total energy density it 
is straightforward to find 



(35) 



When the condensation process ends, for t > IbeCi the evolution of the Universe is given by 



equation (20). As can be seen from the equations (32) and (35) the scale factor do not depend on r 



With equations (28), (32) and (35) we calculate the Hubble parameter at the critical point for the 



model with and without cosmological constant, respectively 



H, 



crt 



(0) 



H, 



crt 



FoVn$ 0) (i + *crt) 3 + nr + nA > 



(0) 



(36) 
(37) 



which allow us to constraint the critical redshift z cr t, since the Hubble parameter at the critical point 
must be bigger than the one today. In Figure ((!]) we see the ratio between the Hubble parameter at 
the beginning of the condensation process and the Hubble parameter today as a fuction of the critical 
redshift z cr t- The requirement that H cr t/Ho > 1 is satisfied for z cr t > 3.5, z cr t > 3.75 and z cr t > 
respectively. 

The effective equation of state parameter u is 



P_ 

pc 2 



a 2 4 0) (l + ^r 3 -Ovi? 3 



of' + ^fi + ^j-H^ 



,(o) 



(38) 



where R = a/ao- This expression and the conservation equation Q allow us to calculate the speed of 
sound c 2 = dp/ dp, that is 

' i \ \ ' . ii — :— / i 

(39) 



C7 2 4 0) (l + ^t)- 3 



C 2 ^ 0) +4 0) (l + Z crt )-3 + O Aj R 3 



which we can approximate to null, since a 1 « 1. 
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Figure 1: The ratio between the Hubble parameter at the beginning of the condensation and the 
Hubble parameter today for a Universe filled with cosmological constant (dashed line) and without 
cosmological constant (solid lines). 



3 Cosmological density perturbations 

There are other studies about this subject: with Newtonian gravity [2D], with the new-Newtonian 
gravity |21|. 122] . with Post-Newtonian approximation by using the conservation of the general rela- 
tivistic energy-momentum tensor, and considering the small velocity limit [23] and the relativistic case 
|24j where the model unifies the dark matter and dark energy. 

Here we study the evolution of the structure formation during the BEC phase transition using the 
gauge-invariant formalism |25[ [26] for the classical linear perturbations. Whe choose the longitudinal 
(or conformal-Newtonian) gauge, in which for a diagonal energy-momentum tensor the perturbed 
metric is 

ds 2 = a(7]) 2 [(1 + 2<S>)d7] 2 - (1 - 2<5>)dx 2 ] , (40) 

where is the gauge-invariant scalar perturbation, which can be interpreted as the relativistic gen- 
eralization of the Newtonian potential and r] is the conformal time, related to the cosmic time t by 
cdt = adrj. We replace the metric (40) in the Einstein field equations ^ and keeping only the first 



order terms we get the equations of motion for the scalar potential <I> 

V 2 $ - m (*' + U<f>) = ^a 2 G5p , (41) 

$" + m& + (m! + v 2 ) <3? = akGcl 2 5p/c 2 , (42) 

where the apostrophe ' means the derivative to the conformal time rj, % = a' /a is the conformal 
Hubble parameter and bp and bp are the gauge- invariant perturbations of the total density and the 
total pressure, respectively. 

For the case of adiabatic perturbation [251 EE] we use the relation 

bp = c 2 s bp , (43) 
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combined with the equations (41), where c 2 = dp/ dp is the fluid sound speed, to find 

r 2 s 



+ m 1 + -4 *' + 



■m' + h 2 {i + + ^k 2 



$ = 



(44) 



where k 2 is the wave-number and we made V 2 $ = — k 2 <&. With help of a change of variable [26J and 
the background equations we find 



,d 2 <S> 
~da 2 



+ 



4 + 3 ^f- \( l +Z—2 
c z 2 \ pc l 



da 



1 + 3- 



1 + 



A' 2 



8-KGpa 2 



$ = 



(45) 



To examine the density perturbations 5p we calculate the relations T = 
; is the covariant derivative, to find |25|, [26] 

V + ^(<5p + 5p) - 3$'(p + p) + a(p + p)^ = , 

a~ 4 [a 5 (p + p)5u\ l ]' + V 2 <5p + (p + p)V 2 $ = , 



2 \ 

up to first order, where 



(46) 
(47) 



where 5u % is the fluid velocity perturbation and , i is the ordinary spatial derivative. 

In this analysis the baryons and DM do not interact except through gravity and we neglect the 



radiation contribution. So each component satisfy the pair of equations (46) and (47) separately. Since 
Pert/ P\C 2 = & 2 « 1 is true we can make an approximation to find that both baryons and DM satisfy 
the same differential equation 



(5-3$)' + a8u\ 
s= Spb = $P X 

Pb P X 







(48) 
(49) 



where 5 is the density contrast. 



Combining equations (47) and (48) and making a change of variables for R = a/ao we find that 

V 2 $ 



d 
dR 



R 2 n^- (5 - 3$) 
dR 



H 







(50) 



From equation ( 38 ) we see that 



n A = 





-n A R 3 



(51) 



where fL 



^+4°)(i+z, 



, , i j " 3 . From equation (39) for the sound speed we have 

'2 



(52) 



since we consider a 2 « 1. 



3.1 Universe filled with baryons and DM 



With equations (74) and (52) the perturbation equation (45) will be 

,d 2 <J> 



R 



dR 2 2 dR 



1 



$ = 



(53) 
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which has a simple power law solution 



n± 



C+R n + +C-R n - 
-5± V33 



(54) 
(55) 



and the growing and decaying modes of the potential $ as a function of time are shown in Figure ^ . 
Applying the solution ( 54 ) to the equation ( 50 ) it is straightforward to find 



3$(i?) + 



tti 



R {n+-2) i?( n -" 2 ) 



n + {n + - 2) 



n_(n_ — 2) 



(56) 



The growing and decaying modes of the matter density contrast as a function of time are shown in 
Figures ^ and Q. 
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Figure 2: Growing and decaying modes of the potential $ as a function of time in a Universe filled 
with baryons and DM. 



3.2 Universe filled with baryons, DM and cosmological constant 

In the case of Universe dominated by a cosmological constant with density parameter Q A we find, 



from equations (45), (51) and (|52j), that 



dR 2 



+ 



1 



1 



n A R 3 \ 



n A R 3 J 



d<S> 



1 - 3 



n A R 3 \ 



n A R 3 J 



and introducing the new variables 



x 



n A R 3 



n m + n A R 3 ' 

x u± y(x) , 



we will have 



d 2 y 

x(l - x )~j^2 + I C ± ~~ + ^± + l) x 



dy 
dx 



(a±b±) y = 



R = 



(57) 



(58) 

(59) 
(60) 

(61) 
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Figure 3: Growing and decaying modes of the matter density contrast as a function of time for various 



wavelengths in a Universe filled with baryons and DM, where Xh = cH 1 is the Hubble wavelength. 
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(62) 

(63) 
(64) 

(65) 



The solution to the equation (57) is 



$(R) = C + x u + 2 F 1 (a + ,b + ;c + ;x) + C-x u ~ 2 F 1 (a_,b_;c-;x) , 



(66) 



where 2^1(0, b; c; x) is the hypergeometric function. The growing and decaying modes of the potential 
$ as a function of time are shown in Figure ([5]). 

In this case we can solve equation ( 50 ) numericaly with help of expression ( 66 ) to find the density 
contrast. The growing and decaying modes of the matter density contrast as a function of time are 
shown in Figure ([6]). 
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Figure 4: Growing and decaying modes of the matter density contrast as a function of time for various 
critical redshifts in a Universe filled with baryons and DM. 
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Figure 5: Growing and decaying modes of the potential $ as a function of time in a Universe filled 
with baryons, DM and cosmological constant. 

4 Cosmological gravitational waves 

In this section we study the gravitational waves generated by the expansion of the Universe during 
the phase transition of normal DM into BEC DM. To do so it is necessary to perturb to first order 
the metric tensor and we write 



ds = c dt — [a(t)5ij + hij(t,x)] dx l dx J 



(67) 



where hij(t,x) are the metric perturbations and we adopt the transverse and traceless (TT) gauge 
[23 [23 EH]) = and ■ = 0. Applying the metric (67) to the Einstein field equations 

Q keeping only the first order terms [25] and with help of the equations of motion ^ we find the 
dynamics equations of the cosmological gravitational waves 



a 



a? a 



hij = 



(68) 
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Figure 6: Growing and decaying modes of the matter density contrast as a function of time for various 
parameters in a Universe filled with baryons, DM and cosmological constant, where \h = cHq 1 is 
the Hubble wavelength. 



where V 2 is the Laplace operator. We introduce the expansion 



h i:j (t,x) = ^ 



d 3 k 
(2tt) 3 /2 ( 



(69) 



where s 



are the two polarization states of the gravitational waves, k is the wave-number and 



elj'(k) is the polarization tensor, where e[f (k)e^ s '^ (k 1 ) = 25 ss /5 3 (k — k') and = = h(t,k). 



Writing the wave equation (68) as a function of the scale factor a, using the (69) and the background 
equations §5§ we find 



da 2 



3 . . dh 



3k 2 



1 



h = 



(70) 



with p and p being the total pressure and density respectively. 



With the solution of the wave equation ( 70 ) we can calculate de gravitational waves energy density, 



that is related to the perturbation hij(t,x) by [27] 



Pgw 



1 



32vrG 



(71) 
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where (• • • ) means spatial average over many reduced wavelengths. The logarithmic energy density is 

_ 8ttG d Pgw 

9W ~m 2 dink ' [{ ) 

and the spectral density Sh(k) [27] is related to the logarithmic energy density by 

n gw = ^k 3 S h (k) . (73) 

4.1 Universe filled with baryons and DM 



In order to find a solution to the wave equation (70) we first consider a Universe filled with baryons, 
which are pressurless, and DM in phase transition. Since a 2 « 1 and, in this case, Q\ = 0, we can 
make the following approximations 

u « , (74) 
rwn x . (75) 



And the wave equation (70) will be written as 



od 2 h 3 dh ( k 2 1\ 

R m-2 R iR + {4BiaZ a+ ^ h=0 



where R = a/a^ and 

° -S, 



n m = n { ° ) + n^(i + Zcrt )- 3 . (77) 



To solve the above differential equation for the gravitational waves it is necessary to introduce the 
variables 

<? 2 = 3^" . (78) 
oo m 

x 2 = ARq 2 , (79) 
h = x 5/2 y(x) , (80) 

and we obtain the Bessel differential equation 

x 2 y"(x) + xy'{x) + (x 2 - 17/4)y(x) = , (81) 

which gives the solution 

h(2q^/R) = (2qVRY /2 (c x J v (2qVR) + C 2 J-v(2qVR)) , (82) 



where J v {x) is the Bessel function of first kind with order u = yl7/2 and C\ and C2 are integration 
constants. In Figure ([7]) we show the logarithmic energy density for several values of the parameters 
for a Universe filled with baryons and CDM during the transition phase as a fuction of the ratio 
between the gravitational waves wavelength and the Hubble wavelength, defined as \h = cHq 1 . 
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Figure 7: Gravitational waves energy density as a fuction of the wavelength for fixed time in the 
case of a Universe filled with baryons and DM during the phase transition, where Xh = cHq 1 is the 
Hubble wavelength. 



4.2 Universe filled with baryons, DM and cosmological constant 



In this section we will consider a Universe filled with baryons, DM in phase transition and a cosmo- 



logical constant. In the equation of state parameter (38) above the first term can be discharged, since 
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a is very small and we can make the approximation 

_ n A R 3 

U} ^~n m + n A R 3 ' (83) 

where 

n m = n° b + n^{i + z crt )~ 3 , (84) 



and the wave equation (70) is written as 

2 d 2 h 3 f n m \dh f q 2 R if n A R 3 \\ h= 

dR 2 2 \n m + n A R 3 )dR \n m + n A R 3 2\ n m + n A R 3 )) ' 1 j 

with 

1 2 = WZ2 ■ (86) 

For the case of perturbations with wavelength bigger than the Hubble length we approximate q ~ 
to find 

2 d 2 h 3 / n m \dh i f n A R 3 \ 
dR 2 2 \n m + n A R 3 J dR 2\ n m + n A R 3 J ' 1 j 

and we introduce the variables 

* = ~R 3 > (88) 
h( x ) = x u± y{x) , (89) 

where v± = 5/12 ± yYJ /12 and we find 

x(l - x) + [c± - (a± + b± + l)x] ^ - a±b±y = , (90) 
dx z dx 

with 



(3-2V5±VT7) 
o± = ^2 ' ^ ' 



b± _ (3^ (92) 



V17 , x 

c± = 1±^- • (93) 

The solution is found to be 

/i(i?) = C + x y + 2 F 1 {a + ,b + -c + -x) + C-x u - 2 F 1 (a-,b^c^x) , (94) 

where 2-Fi( a 5 b; c; x) is the hypergeometric function. 
For small wavelenghts, the equation (85) will be 

2 d 2 h 3 u f n m \dh q 2 R 



2 \n m + n A R 3 JdR n m + Q A R 3 

which can only be solved numericaly. In Figure Q we show the spectral density, that in this case 
does not depend on the wavelength, as a function of time in a Universe filled with baryons CDM and 
cosmological constant during the transition. 



16 



ti b 0) = 0.045, n <0) =0.255, n A =0.7 (C 2 =0) 



8x 10" - 




ti b 0) = 0.045, n <0, =0.255, n A =0.7 (Ci=0) 



8x 10 16 - 



0.000 0.001 0.002 0.003 0.004 0.005 

t 

to. 




Figure 8: Gravitational waves spectral density as a fuction of time in the case of a Universe filled with 
baryons and DM during the phase transition. The black curve shows the case where z cr t = 7 and in 
the gray one z cr t = 12. 

5 Conclusions 



In this letter we considered a Universe where initially the DM is nearly pressureless and after a phase 
transition the DM is all in form of BEC DM, which has a polytropic equation of state. During 
the transition phase both non-condensated and condensated DM coexist and the time evolution of 
the Universe expansion is changed. We studied the scalar and tensorial perturbations generated by 
the background expansion during the phase transition in two cases: with and without cosmological 
constant. 

Firstly, we verify from Figure [l] and from equations (36 37) that the ratio between the Hubble 
parameter at the beginning of the condensation process and the Hubble parameter today, with the 
condition H cr t/Ho > 1, is satisfied for z cr t > 3.5, z cr t > 3.75 and z cr t > 0, what was expected. 
Two situation are analysed: for a Universe filled with cosmological constant and without cosmological 
constant. We can see that the case where A is present is favored in the sense that the ratio of the 
Hubble parameters is bigger than the case where the cosmologcial constant is absent. This indicates 
that in a Universe whose material content is baryonic matter, BEC DM and cosmological constant 
the condensation can begin later when compared with the model without A. 

When we analyze the scalar potential $ we see that in the case without cosmological constant, 
Figure [2j there is a growing and a decaying mode. Both modes are more sensitivy to changes in the 
amount of DM than to changes in the critical redshift z cr t- In the Universe with cosmological constant, 
Figure [5j we see two decaying modes and smaller the critical redshfit z cr t smaller the potential 

In Figures [3] and Figure [4] we see that the density contrast amplitude, in a Universe without 
cosmological constant, is sensitivy to the amount of DM and bigger the perturbation wavelength 
bigger the perturbation in the case of the growing modes and smaller in the case of the decaying 
modes. The same behavior can be seen when the critical redshfit is varied. 

For a Universe with cosmological constant the contrast amplitude in Figure [6] can not be distin- 
guished when we vary the perturbations wavelength in the case of the decaying modes. For the growing 
modes the amplitude initially decays but after some time interval it grows again. The decaying period 
of the growing modes depends on the wavelength. The bigger the wavelength bigger is the decaying 
time interval. The same behavior is reproduced when W6 vary trie critical redshift z cr f. 
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In the case of tensorial perturbations in a Universe without cosmological constant we see in Figure 
[7] that the gravitational waves energy density changes little when we vary the time of beginning of 
the transition (z cr t) and even the amount of DM in the Universe. Figure [8] shows little changes in 
the spectral density in a Universe with cosmological constant when we change the beginning of the 
transition. 

The Bose-Einstein condensated dark matter model can be an interesting theory to explain the 
nature of dark matter in the Universe. A more detailed investigation of the model discussed herein, 
including the generalization of the equation of state of BEC dark mater model, is being made as well 
as the bayesian analysis to determine the most probable values of the various parameters of the theory. 
The BEC DM model still has to pass by several observational tests to prove as a viable cosmological 
model. Its importance grows as weakly interacting massive particles, or WIMPs, are searched in the 
Large Hadron Collider or other ground-based experiments. 
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